Driving light pulses with light in two-level media 
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A two-level medium, described by the Maxwell-Bloch (MB) system, is engraved by establishing a 
standing cavity wave with a linearly polarized electromagnetic field that drives the medium on both 
ends. A light pulse, polarized along the other direction, then scatters the medium and couples to 
the cavity standing wave by means of the population inversion density variations. We demonstrate 
that control of the applied amplitudes of the grating field allows to stop the light pulse and to 
make it move backward (eventually to drive it freely). A simplified limit model of the MB system 
with variable boundary driving is obtained as a discrete nonlinear Schrodinger equation with tunable 
external potential. It reproduces qualitatively the dynamics of the driven light pulse. 

PACS numbers: 42.50.Gy, 42.65.Re 



Introduction. Manipulation of light with light has be- 
come one of the hottest research spots in quantum op- 
tics this last decade. A widely studied field of research 
makes use of electromagnetically induced transparency 
in three-level systems, which allows to slow down, and 
eventually stop, a light pulse [1, 2, 3]. Another intersting 
research option uses resonantly absorbing Bragg reflectors 
(RABR) which consist in a periodic array of dielectric 
films separated by layers of a two-level medium A, 5,, Q , 
allowing a light pulse not only to be stopped and trapped 
01 , but also to be released by scattering with another 
control pulse, thus creating a "gap soliton memory" Q. 

The fundamental process underlying such novel light 
pulse dynamics is the cooperative action of nonlinearity 
and periodicity. Still, a serious drawback when making 
use of RABR is the built-in periodic structure that re- 
stricts both the freedom of pulse parameter and of pulse 
dynamics. 

We propose to prepare a two-level system (TLS) by 
establishing a standing electromagnetic wave in a given 
polarization direction, and then to scatter a light pulse, 
orthogonally polarized. The incident pulse then feels the 
electromagnetic induced grating through the coupling me- 
diated by the population density, as described by the gov- 
erning Maxwell-Bloch (MB) system [i,!!^. The freedom 
in the choice of the standing wave parameters (in partic- 
ular the boundary amplitudes) allows us to demonstrate 
by numerical simulations as in Fig[T] that the incident 
light pulse not only can be stopped but also can be re- 
leased back to the incoming end. Engraving a medium 
with a cavity standing wave is a method previously used 
to create two dimensional w aveg uide arrays in strongly 
anisotropic photonic crystals 

In a TLS of transition frequency il, the MB system 
is considered in the isotropic case for a plane polarized 
electromagnetic field propagating in direction z. The 
time is scaled to the inverse transition frequency f^"^, the 
space z to the length U,c/rj {jj is the optical index of the 
medium), the population inversion to the density of ac- 



tive dipoles A^o, the energy to the average Wq — Nohn/2 , 
the electric field to y^Wo/e and the polarization to \/eWo. 
The resulting dimensionless MB system then reads 



Ett — -1- Pu 
Ptt+P + aAE = -72Pt, 
Nt-E-Pt^ -71(1 + A). 



(1) 



where E and P, denote vectors in the tranverse plane, 
e.g. E — {Ex{z,t), Ey{z,t)). The coupling eventually 
results in a unique dimensionless fundamental constant 
a = 2fioc^/{3r]^)\ni2\^No/ihn){r]^ + 2)^/9 (where the 
dipole moment /ii2 is averaged over the orientations [lol |) 
and by the normalization of the population inversion den- 
sity: A^ = — 1 when all active atoms are in the fundamen- 
tal state, A^ = 1 in the excited state. The dimension- 
less dissipation coefficients are 7 — {c/r]){A/^) resulting 
from the electric field attenuation A, then 71 = l/(f7Ti) 
and 72 = 2/(^1X2) resulting respectively from the popu- 
lation inversion dephasing time Ti and the polarization 
dephasing time T2. 

In the strong coupling case a '--^ 1 (dense media), a 
multiscale analysis has shown that the model equation for 
the two directions of polarization results to be a system 
of coupled nonlinear Schrodinger equations Then 
reduced to a unique polarization, it allows for slow light 
soliton formation out of evanescent incident light [l3|, 
namely under irradiation in the forbidden bang gap. This 
gap results from the linear dispersion relation of the MB 
system, for a carrier exp[i(wt — kz)] on a medium at rest, 
namely 
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The upper edge of the stop gap at frequency ujq corre- 
sponds to fc = and du/dk = 0. 

Numerical simulations. We first proceed with numer- 
ical simulations of the MB equations ^ submitted to the 
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following boundary-value problem 



E(0,i) 



E(L,<) 



£2it) C0s{LLj2t) 
f^'^^ COs((jJit) 





(3) 
(4) 



where £2 (t) is the slowly- varying low-amplitude pulse en- 
velope 
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cosh[/i(f — fo)] 



(5) 



The carrier frequencies are chosen close to the gap edge 
wqj inside the passing band for the grating field, in the 
gap for the incident pulse, namely 

wi > Wo, < t^o, - t^ol ^ C'(e^) (6) 

for J = 1,2, where e is our small control parameter. We 
set from now on a = 1. 
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FIG. 1: (Color in line) A typical result of a numerical simu- 
lation of IT]) under boundary values (|3]|4]) with parameters 

^(0) ^ ^(L) ^ ^ ^™a. ^ QQ^^ ^ ^3^^ 

to = 1000, fi — 1/150. The plot shows the evolution of the 
energy density along the polarization direction y, trapped by 
the grating induced by the stationary wave. The pulse is even- 
tually released back at time 4500 by increasing the amplitude 
on the incident side by 25%. 

We have chosen damping coefficients such as to reach a 
fast stabilization of the lattice created by the stationary 
driving: in the x-direction we have set 7 = 72 = 0.01 
and 7i — 0.001, while in the y-direction they have been 
set to zero in order to see clearly the process. Under 
such boundary data, with the parameter values indicated 
in the caption, we obtain the result displayed in FiglT] 
More precisely, the applied grating field (polarized along 
x) is settled smoothly and, at time t = 1000 the incident 
pulse enters the medium. It is naturally trapped and seen 



to oscillate about the center z = 30. At time t = 4500 
the left-hand-side amplitude fj*'"' is increased from 0.4 
to 0.5, which produces the reverse motion of the stored 
pulse. Actualy one may paly with the driving amplitudes 
and £[^^ to drive the pulse back and forth. 
Interpretation. The process described above is now 
understood, within a multiscale analysis of MB equa- 
tions, in terms of a discrete nonlinear Schrodiger model 
with variable coefficients related to the variation of the 
boundary grating field amplitudes. Thanks to (O, a so- 
lution of ll]) under boundary values ([5]|3]) can be sought 
under the form 



E(z,t) 



e^E2{C,Ti,T2) 



(7) 



with the slow variables ( — ez and t„ = e^"t. The second 
slow time T2 is meant to capture the nonlinear dynamics 
of the low-amplitude (e^) long duration (e^^) incident 
pulse. Note that form assumption ([6]), the frequency shift 
from 0^0 is contained in the slow time variations with 
n = e^t. 

Inserting then fT]) in the MB system ([T]) we eventually 
obtain (with vanishing damping) 
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where star means complex conjugation. The boundary 
values ([3114]) for the x-component Ei imply from equa- 
tion ([8]) that £'i(C, Ti) is a periodic stationary solution 
of the nonlinear Schrodinger equation with frequency 
vi = {uji — ciJo)e^^- It acts then as an external periodic 
potential in the evolution 

In order to take into account variations |i?i|(C) result- 
ing from the variations of the boundary driving ^ we 
set (remember viTi = (wi — 0Ja)t) 



i?i-|i?i|(C) e'' 



m^O^V.iO+e'ViC) (10) 



where Vo(C) is purely periodic while V{C,) describes 
the aperiodic inhomegenities of |i?i|(C) induced by the 
boundary values. 

We then seek a solution of ([5]) on a suitable orthonor- 
mal basis of Wannier functions Lpj {Q) which are localized 
with respect to the site j [l3|, within the one-hand ap- 
proximation [ist , where j actually indexes the minima of 
the periodic potential Vo(C)- A solution of ^ is sought 
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^j(n,T2,)(/3j(C), 



(11) 
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and the equation for the coefRcients is worked out by 
inserting (jlip in ([9]) and by projecting on a chosen (pj. In 
the tight-binding approximation [l6| we eventually obtain 



- k,T] + Q(^,_i + + C/|^,f , (12) 

where the coefficients are given from the Wannier basis 
by (integrals run on C € M) 
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Note that translational invariance guarantees that the 
above coefficients are j-indepcndcnt, except of course for 
U,j and Aj that bear the aperiodic inhomogcncity of the 
external potential 

Equation (|12p can now be solved first for the ti- 
dependence of Tj as a linear system, which provides then 
a discrete nonlinear Schrodinger coupled system for the 
T2-dependent amplitudes. This is done by seeking a so- 
lution under the form 



1 lAri 



(13) 



for which the leading order of furnishes the two cou- 
pled linear equations 
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The dispersion relation and the relation between Qj and 
(^G^^)* automatically follows as 



Ao 



1^1 



At next order we readily get for Qj the discrete non- 
linear Schrodinger equation, the parameters of which 
are greatly simplified if we note that in our numerical 
simulations (for sufficiently deep lattice) f2o 3> Aq and 
ilj ^ Aj . In such a case the equation reads 



n,Gj - Q[Gj+, + + U\Gj\^Gj. (14) 



The electric field envelope in the y-direction of polariza- 
tion reads 



E 



(15) 



in terms of the solution of the chosen Wannier basis. 

The dynamics of the pulse is thus interpreted out of the 
discrete nonlinear Schrodinger model (fT4|) as the action 
of the potential Vlj that translates the applied variations 
of the boundary driving in the cc-direction of polarization. 

Application. In order to illustrate the above inter- 
pretation, we proceed now with numerical simulations 
of ([M]) where T2 = t to read quantities in physical di- 
mensions. The potential Qj models the variations of 
\Ei\ away from a purely periodic function as soon as 
^£[^\ We set 





'0.07 J 



t < 650, 
t > 650 



(16) 



and obtain the Figl2] which shows the same qualitative 
behavior as Fig[T] 
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FIG. 2: (Color in line) Plot of the solution of HHl with the 
parameters Q = 0.2, U = 0.4 and an incident pulse Go{t) = 
0.305e"° ''^*"/cosh[(f-750)/250] which propagates, becomes 
trapped and is then released back by the action of flj in (fT6]) . 

Conclusion and comments. We have shown that the 
Maxwell-Bloch system, the celebrated fundamental and 
general semi-classical model of interaction of radiation 
with matter, may serve as a tool to store and drive light 
pulses in an arbitrary way by conveniently using laser 
light to engrave the medium with a controllable stand- 
ing wave pattern. The extreme genericity of MB model 
together with the freedom in the boundary values of the 
engraving field, constitute a decisive advantage over other 
techniques to store and manipulate light pulses. 

The process is understood by deriving a discrete non- 
linear Schrodinger model where the external tunable po- 
tential actually translates the effects of boundary driving 
variations. Our purpose was simply to provide a qualita- 
tive interpretation. 

A more detailed study, reported to future work, would 
require first to use exact solutions of ([H]) for the grating 
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field (e.g. in terms of Jacobi elliptic functions), second 
to construct the corresponding most adequate Wannier 
basis, third to evalute precisely the effect of boundary 
driving variations on the factor f2j , and last to study the 
full system that couples to T*. 
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